Eigenvalue asymptotics of perturbed periodic 
Dirac systems in the slow-decay limit. 



Karl Michael Schmidt 
School of Mathematics, Cardiff University, 23 Senghennydd Rd 
Cardiff CF24 4YH, UK 

Abstract. A perturbation decaying to at oo and not too irregular at introduces 
at most a discrete set of eigenvalues into the spectral gaps of a one-dimensional Dirac 
operator on the half-line. We show that the number of these eigenvalues in a compact 
subset of a gap in the essential spectrum is given by a quasi-semiclassical asymptotic 
formula in the slow-decay limit, which for power-decaying perturbations is equivalent 
to the large-coupling limit. This asymptotic behaviour elucidates the origin of the dense 
point spectrum observed in spherically symmetric, radially periodic three-dimensional 
Dirac operators. 

1 Introduction. 



For a large class of potentials, the semiclassical Weyl formula gives a correct asymptotic 
description of the total multiplicity of the lower spectrum of a Schrodinger operator 
in the large coupling limit (see [18] Chapter XIII. 15 and the references given there). 
However, there are also some notable exceptions, e.g. in the two-dimensional case ([9], 
[10]). In recent years an analogous asymptotic analysis of the point spectrum arising in 
spectral gaps of Schrodinger operators under perturbations has attracted considerable 
attention; starting from [1], Birman has developed a general framework to study this 
problem [2], [3], [4], [5], [6], [7]; see also [16]. 

More specifically, Sobolev [26] has studied the perturbed periodic one-dimensional Schro- 
dinger operator, showing that for a wide range of power-decaying perturbations, the 
number of eigenvalues arising in a closed subinterval of a spectral gap of the unper- 
turbed problem is asymptotically given by a quasi-semiclassical formula in which the 
quasimomentum of the periodic background problem takes the role of the ordinary mo- 
mentum in the Weyl formula. This result has an interesting application to spherically 
symmetric Schrodinger operators in BrC with a radially periodic potential, which have 
dense point spectrum in all spectral gaps of the corresponding one-dimensional peri- 
odic operator [14]. Indeed, treating the angular momentum term of the one-dimensional 
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half-line operators in the partial-wave decomposition 



-A + g(|.|)=© -^+q{r) + {l{l + n-2)+'-^^-^f—^)/r^ (r G (0, oo)), 



with periodic q, as a perturbation, the quasi-semiclassical formula explains the origin 
of the dense eigenvalues, and describes their asymptotic density in the limit I — > oo. 
Furthermore, numerical experiments have shown that the asymptotic formula can give 
a surprisingly accurate prediction for the number of eigenvalues even for small values 
of the coupling constant [11]. (We mention in passing that the angular momentum 
term also happens to be a boundary case for the question whether the total number 
of eigenvalues introduced in a spectral gap for an individual half-line operator is finite 
or infinite ([19], [20], [21], [22], [25]); but this critical behaviour is not apparent in the 
large-coupling limit.) 

The relativistic counterpart of the Schrodinger operator, the Dirac operator, is un- 
bounded below, so that there is no lower spectrum, and one is always in a gap situation 
when studying the discrete spectrum, even when there is no periodic background poten- 
tial (for the large coupling asymptotics in that case, see e.g. [8], who correct a result by 
[15], and [12]). Again, an interesting class of perturbed periodic one-dimensional oper- 
ators arises from the partial-wave decomposition of the spherically symmetric, radially 
periodic Dirac operator in IR^, 

d k 
H = -ia ■ V + ml3 + q{\ ■ \) ^ ffi -ia2— + mas + q{r) + -cri 

dr r 

where ai, a2, cx^, (3 are 4x4 Dirac matrices, and 



This operator also has dense point spectrum in the spectral gaps of the corresponding 
one-dimensional periodic operator, but one exceptional gap may or may not contain 
dense point spectrum — it is an open question whether partial filling of the gap can 
occur [23], [24]. The underlying reason for the complications in the Dirac case is the 
fact that the angular momentum term aik/r, in contrast to the Schrodinger angular 
momentum term, docs not have a definite sign. Since it does not depend monotonically 
on the quantum number k, the eigenvalues can move either way as k increases. Methods 
such as Sturm comparison, which yields a fairly quick proof for the quasi-semiclassical 
formula for perturbed Sturm-Liouville operators, are of little use when applied to matrix 
perturbations of Dirac operators. 

In the present paper we study the asymptotic distribution of eigenvalues in spectral 
gaps of periodic Dirac operators with a matrix perturbation of the type of the angular 
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momentum term. We assume that the perturbation tends to at oo, but no assumptions 
on its sign or rate of decay are needed. The asymptotic Umit we consider is in fact a 
'slow-decay' Umit, which for perturbations of inverse power decay is equivalent to the 
'large-coupling' limit, an observation central to Sobolev's result. More explicitly, given 
a Dirac system 

■ d 

T = + maz + q[x) 

with m > and real- valued, a-periodic q G Ll^^{IR), and a real- valued perturbation 
template 

Zo £ C{0, oo), lim Zo(^) = oo, lim Io{q) = 0, 
we investigate the number of eigenvalues of the perturbed Dirac system on (0, oo) 

r(c) = -i^2^ + maz + q{r) + lo{r/c)ai 

in a compact subintcrval of a spectral gap asymptotically as c ^ oo. 
The angular momentum term is clearly of this type, with Io{q) = 1/q. 
We follow Sobolev's basic idea of comparing t(c) to an operator in which lQ{r/c) is 
locally replaced by a constant, but in the absence of Sturm comparison to obtain upper 
and lower bounds we must resort to the coarser instrument of operator perturbation 
theory. This complication leads to severe difficulties near 0, where the perturbation is 
divergent: in the Sturm-Liouville case this actually helps, as the perturbed eigenvalue 
equation becomes disconjugate near 0, so that no spectrum is produced near that end- 
point. The case of the Dirac system, however, with an unperturbed operator unbounded 
below and a strong perturbation of no definite sign, is entirely different. For this reason 
we introduce the following additional assumptions on q and Iq to ensure that essentially 
no spectrum is created near 0; as will be seen, they also imply that r(c) is in the limit 
point case at for sufficiently large c. 



{H) q e L°°{my, lo e ACiociO, lim sup < oo. 



\l'oi6)\ 
1-^0 llio) 



The requirements on Iq are clearly satisfied in the case of the angular momentum term. 
Let k{\, I) be the quasi-momentum (cf. Section 2 below) of the periodic equation 



Tu = {\ — lai)u, 



for x,i e m. 

Our main result is 



Theorem 1. Assume that (H) holds in addition to the general hypotheses, and let 
[Ai, A2] be compactly contained in a spectral gap of the self-adjoint realisation h of t. 
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Then for sufficiently large c > 0, t(c) is essentially self-adjoint on Cq°(0, oo), and the 
number N{c) of eigenvalues in [Ai, A2] of its self-adjoint extension h{c) satisfies 

N(c) 1 f°° 
lim ^ = — / ik{X2, loiQ)) - HXu loiQ))) dg. 
c—>-oo c air Jq 



Remarks. 

1. Our assumptions imply that the essential spectra of h{c) and h coincide, so /i(c) 
only has discrete eigenvalues in [Ai, A2]. It will be apparent from the proof that the 
asymptotic formula in Theorem 1 continues to hold if, instead of (H) and lim lo{g) = 00, 

we assume £ ^"([0, 00)). 

2. In [23] Theorem 2, it was shown that A G iR is a point of the essential spectrum of H 
if A is in the spectrum of the periodic operator h + lai for some / G IR, i.e. if A is a point 
of growth of fc( • , i). In a rather weak sense. Theorem 1 is a reverse of this statement: if 
a A interval does not intersect the spectrum of ^ + lai for any I G IR, then the integral 
on the r.h.s. of the asymptotic formula vanishes, i.e. the asymptotic eigenvalue density 
is 0. Of course, this does not rule out the existence of dense point spectrum of H in this 
interval. 



2 Rotation number and quasimomentum of periodic Dirac systems. 



The qualitative behaviour of solutions of a linear ordinary differential equation with 
periodic coefficients is characterised by its monodromy matrix, the value of the canonical 
fundamental system after one period. In the case of an a-periodic Sturm-Liouvillc or 
Dirac equation, this is a 2 x 2 matrix with real entries and (Wronskian) determinant 
1; the position of its eigenvalues ^1, 1x2 in the complex plane is thus fully determined 
by its trace, the discriminant D = jii + //2- If / ^2, then there are corresponding 
Floquet solutions -ui, U2 satisfying 

Uj{x + a) = jij Uj{x) {x e M, j E {1, 2}). 

Hence it is easy to sec that the equation is stable (all solutions arc globally bounded) 
if |D| < 2, and unstable (all solutions are unbounded) if > 2; as D is an analytic 
function of the spectral parameter, the real line splits into alternating stability and insta- 
bility intervals. This already provides a complete description of the (purely absolutely 
continuous) spectrum of the corresponding self-adjoint ordinary differential operator: 
it is the closure of the union of all stability intervals, spectral gaps corresponding to 
non-degenerate instability intervals (cf. [13], [27]). 
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Furthermore, the oscillation behaviour of solutions is very closely linked to the discrim- 
inant as well. Thus for the periodic Dirac system, it is not difficult to verify along the 
lines of [13] proof of Thm 3.1.2, that for any value of the spectral parameter in the 
closure of the nth instability interval, n G the Priifer angle '& of any iR^-valued 
solution u, defined by 

u = R 



with R > 0, satisfies 



cos ■& 
- sin 



Mx) = !^ + o(l) (x ^ oo). 
a 



In particular, the asymptotic rate of growth of i?, the so-called rotation number 

iim ^^-^ = — , 

is constant in instability intervals. 

It turns out that the rotation number is in fact well-defined as a continuous non- 
decreasing function of the spectral parameter on the whole real line, which can be 
expressed in terms of the discriminant in the stability intervals. For periodic Sturm- 
Liouvillc equations this was elegantly shown in [17], using a connection between the 
winding number of a complex-valued Floquet solution in the punctured plane and the 
Priifer angle of its real part. Unfortunately, this argument does not carry over to the 
Dirac system, with Floquet solutions (in the case of stability) consisting of two complex- 
valued components with no immediate link between them. Nevertheless, the following 
statement holds. 



Theorem 2. Let •& be the Priifer angle of an IR^ -valued solution of the Dirac system 

(-io'2-7-- + rna3 + lai + q)u = Xu 

ax 

with real-valued, a-periodic coefficients m, I, q & Ll^^{IR), X & IR in a stability interval, 
and let D be the discriminant. 

Then there is k{X) G M such that i9(x) = ^^^^^ + 0(1) (x — > oo), and D = 2 cos k{X). 

k{X) is called the quasimomentum of this periodic equation; k{X)/a is the rotation 
number. k{X)/{a7r) is called the integrated density of states in view of the following 
consequence of Theorem 2, which is readily obtained using oscillation theory ([27] Thm 
14.7, 14.8). 



Corollary 1. In the situation of Theorem 2, let {an)neiN, {bn)neiN be sequences 
of real numbers such that an < bn and lim 6„ — a„ = oo, and tn any self-adjoint 
realisation of 

■ d 

—io-2-j h maa -\- lai -\- q 
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on [ttn, bn] with separated boundary conditions. Then for Ai < A2 the number Nn of 
eigenvalues of tn in [Ai, A2] satisfies 

liin ^" = fc(A2) -fc(Ai) 

n^oo hn — an CXTT 



In the proof of Theorem 2 we use the following elementary observation on complex 
solutions of general Dirac systems with real-valued coefficients. 

Lemma 1. Let u : I ^ be a non-trivial solution of 

i-io'2-i- + rnuz + lai + q)u = Xu 
ax 

(m, I, q E -^ioc(-f) real-valued, A G M). Then the following statements are equivalent: 

a) There is x E I, (p E C, \(p\ = 1, such that (pu{x) G . 

b) There is if E G , \ip\ = 1, such that (pu : I ^ IR^ . 

c) u,u are linearly dependent. 



Proof of Theorem 2. 

As \D\ < 2, the monodromy matrix has complex eigenvalues /i 7^ /Z, = 1. By Lemma 
1, the corresponding Floquet solutions u, u are linearly independent as, jj, ^ IR. 
Again by Lemma 1, the components ui, U2 have no zeros, and their arguments (locally 
absolutely continuous functions) are nowhere equal mod tt. By choosing an appropriate 
branch of the argument, we can assume < |arg ui — arg U2\ < throughout. More 
specifically, there are locally absolutely continuous functions Ri, R2 > 0, (fi, {p2, and 
V € 1} such that Ui = i?ie*'^% U2 = ivR2e'^'^^ , and [991 — 992I < ^/2. Interchanging 

u and u if necessary, we can assume without loss of generality that u = 1. 
Consider the iZ?^ -valued solution v := Re u; its Prxifer angle satisfies 

V2 R2 sin f2 R2 , . f. . . , 

tan i} = = — = — cos((/?2 - Vi) (tan ipi -\- tan((^2 - Vi))- 

vi Ri cos (fi Ri 

Thus 'd and (pi arc connected by a modified Kepler transformation (cf. [25]); as a con- 
sequence, they take the values + l/2)7r at the same points, and their difference is 
globally bounded. 

Prom ui{x + q) = /xui(x) we find that i?i is a-periodic, and that there is k E IR such 
that ft = e**^ and (pi{x -\- a) = pi{x) + (x G IR). Thus (pi{x) — is a-periodic and 
continuous, hence globally bounded, and it follows that 

^x) = — + 0(1) {x ^ 00). 
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Furthermore, D = ^ + jj, = 2 cos k. 



□ 



3 Proof of Theorem 1. 

We shall use the following consequences of the spectral theorem (cf. [24] Lemma 6). 

Lemma 2. Let L be a self-adjoint operator with purely discrete spectrum of finite 
total multiplicity Nj[L) in the real interval I = [Ai, A2]. 

(Decomposition Principle) If L is a k-dimensional extension of a restriction of a 
self-adjoint operator L' , then Ni{L) < Ni{L') + k. 

6J (Bounded Perturbations) If A is symmetric and bounded, then Ni{L) < Ni>{L + A), 
where I' = [Ai - A2 + \\A\\]. 

The proof of Theorem 1 proceeds as follows. We introduce the rescaled independent 
variable Q = r/c^ (0, 00), and use the decomposition principle to split off the two regions 
near the singular end-points and 00. The regular operator on the intermediate interval 
is then, by another application of the decomposition principle, split into operators on a 
large but finite number of subintervals on which the perturbation (TiZo(^) does not change 
very much. Lemma 2 b) relates the number of eigenvalues of each subinterval operator 
to that of an operator with periodic coefficients on the same subinterval. By virtue of 
Corollary 1, the asymptotic number of eigenvalues for the latter can be expressed by 
means of the quasimomentum, considering that in terms of the original variable r, the 
subintervals grow beyond all bounds in the limit c 00. 

The operator on the unbounded interval near 00, on which the perturbation is small, is 
easily dealt with in a similar manner. The end-point 0, however, presents more serious 
difficulties, which reflect the lack of monotonicity of spectral behaviour with respect to 
a matrix perturbation; and it is at this point only that our rather strong hypotheses (H) 
enter. We show that with a suitably chosen boundary condition, the operator on the 
interval near has no eigenvalues at all in [Ai , A2] for sufficiently large c. The asymptotic 
formula then follows in the limit of infinite refinement of the subintervals. 

The Dirac system 

r(c)n = \u 

is in the limit point case at infinity ([27] Cor. to Theorem 6.8). We begin the proof of 
Theorem 1 by showing that for sufficiently large c, it is in the limit point case at as 
well, thus proving the essential self-adjointness. 
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Let ^ > and C > be such that Io{q) > {q E (0, g)), and 

<C (.e(0,,,). 

Then (— log loYig) < CIq{q), and hence 

Io{q) < Io{q) exp (^C j'^ loj {q e (0, Q)); 

and similarly 



k{Q) - hq) =- r io<c r ii 

J Q J Q 

which implies /q ^ -^^(0, • ) in view of lim Iq{q) = oo. 

As 771(73 + g — A is bounded, it is sufficient to establish that 

(-«cr2-7- + (^ilo{r / c))u = 
ar 

is in the limit point case at 0. If u is a non-trivial iR^-valued solution of this equation, 
then u'i{r) = —lo{r / c)ui{r) , and setting v{q) := ■ui(c^), we find for c > C 



v{q) = v{q) exp ic J loj > J^^o{q), 



sov ^ L2(0, • ). 

To prove the asymptotic formula, we first observe that the region near does not 
contribute to the eigenvalue count. Indeed, let cq := C + 1; then 

foiQ)--\l'oiQ)\>^^^oo ie^O), 

Co Co 

and thus there is £io > such that 

^o(^) - -KiQ)\ > (Iklloo + max{|Ai|, IA2I} + 1)'. {g G (0, ^o))- 
Co 

Let c > Co, and /io(c) the self-adjoint realisation of r(c) on (0, c^?o) with the boundary 
condition 

UlicQo) + U2{cQo) = 0. 

As r(c) is in the limit point case at 0, /io(c) is essentially self-adjoint on Dq := {u G 
D{ho{c)) I n = near 0}. Hence if A G a{ho{c)) and s G (0, 1), there is u G Dq \ {0} 
such that \\{ho{c) — X)u\\ < £\\u\\, so 

\\-ia2u' + {ma3 + /o( • /c)o"i)n|| < (||g||oo + |A| + 
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But on the other hand, integration by parts yields 

\\-ia2u' + {mas + ^o( • /c)<^i)u\\'^ 

= / (|«f + (m^ + /o( • /cf)\uf - miiaiufu)' + lo{ ■ /c){{asufuy) 
Jo 

= m\u{ceo)? + / (|nf + (m^ + /o( • /cf)\u\^ - ■ /c){{asufu)) 
Jo c 
/•ceo 1 

> / Mr/cf - -\l'o{r/c)\)\u{r)f dr > {\\q\\^ + max{|Ai|, IA2I} + lf\\u 



Thus |A| > max{|Ai|, IA2I}; in particular, h{c) has no spectrum in [Ai, A2]. 

A similar calculation shows that for all fixed Q & {0, qq), the self-adjoint periodic operator 



-^o■2— + m(73 + q{x) + lo{Q)cri 

has a spectral gap, and hence instability interval, containing [Ai, A2], and it follows that 
k{Xi, lo{g)) = k{X2, lo{g)) (g e (0, ^o))- 



We now turn to the remaining interval {ego, 00). Choose S^o G (0, 1) so small that 
Soo < (A2 — Ai)/2, and [Ai — 5oo, A2 + 5oo] is still compactly contained in the spectral 
gap of the periodic problem. Fix Po > go such that |/o(f?)j < 600 {g > ^o)- 
Let e > 0. There is L > such that |/o(^?)| ^ L (g ^ (go, 00)). As the quasimomentum 
k{X, I) is uniformly continuous on the compact set K := [Ai — 1, A2 + 1] x [— L, L], 
there is 5 G (0, S^o) such that 

in, I), in', I') eK,\n- /x'l, \l-l'\<6^ \k{n, I) - k{fj,', < £. 

As Zo is uniformly continuous on [g^ , Pq] , there is 7 > such that 

|/o(a;) - /o(y)| < 5 (a;, y e [go, Po], \x - y\ < 7). 

Now consider a partitioning of the interval / := {go, Pq) into n subintervals Ij = 
{gj-i, gj) with \Ij\ < 7 (j G {1, . . . , n}). 

Let hj{c), J G {1, . . . , n}, and /ioo(c) be self-adjoint realisations of r(c) on clj, and 
(cPq) 00), respectively, with separated boundary conditions at the regular end-points. 

n 

Then ho{c) © hj{c) © hoo{c) is a 2(n -|- l)-dimensional extension of the 2(n -|- 1)- 
i=i 

dimensional restriction of h{c) with domain {u G D{h{c)) \ u{cgj) = 0{j G {0, . . . , n})}; 
comparing these operators by means of the decomposition principle (Lemma 2 a) we 
find 

n 

\No{c) + J2 ^i(c) + ^00 (c) - N{c)\ < 2(n + 1) 
i=i 

for the total spectral multiplicities in [Ai, A2]. As observed above, iVo(c) = 0. 
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To estimate Nj{c) for j G {1, . . . , n}, choose Qj G Ij, let Ij := lo{Qj), and define self- 
adjoint operators 



on c/j- with D{hj{c)) = D{hj{c)). Then Lemma 2 b) and |/o(£') — lj\ < ^ (q ^ Ij) yield 
the bounds 

Nj{\i + (5, A2 - (5; c) < iVj(c) < iVj(Ai - (5, A2 + (5; c), 

where Nj{ni, fi2', c) is the number of eigenvalues of hj{c) in [/xi, 
As /ij(c) has a-periodic coefficients, Corollary 1 implies 



lim 



c—^oo c\Ij\ air 
Thus we find 

^(A;(A2, Ij) - HXi, Ij) - 2e) < ^(fc(A2 - 5, Ij) - k{Xi + S, Ij)) 
air air 

= lim -Nj{Xi + (5, A2 - 5; c) < lim inf 

c— >oo C c— >oo C 

N-(c) 1 ~ 

< lim sup < hm -Nj{Xi - S, X2 + 6; c) 

c— >oo C c— ^00 c 

= ^(A:(A2 + 6, Ij) - k{Xi - 6, Ij)) < ^(A:(A2, Ij) - k{Xi, Ij) + 2e). 
air air 

In order to estimate iVoo(c), we observe that 

Noo{Xl + Soo, X2 - Soo; c) < iVoo(c) < iVoo(Al - Soo, A2 + (5oo; c), 

where iVj (^i, ^^2; c) is the number of eigenvalues in [/ii, /i2] of the self-adjoint realisation 
hoo{c) of r on (cPq, 00) with D{hoo{c)) = D{hoo{c)). Since /loo(c) is a direct summand 
of a 2-dimensional extension of a restriction of the unperturbed periodic operator on 
IR, we find, again by the decomposition principle, that < iV<x)(/ii, 1^2] c) < 2 for any 
[1^1, A*2] compactly contained in a spectral gap. 
Summing up, we obtain 



n n 

— ^(fc(A2, Ij) - MAi, lj))\Ij\ 1/| < lim inf - ^ Nj{c) 

an ^-^ air c ^-^ 

= lim inf ^^'^ < Hm sup = fim sup - ^ Nj (c) 

1 " 2f 

< — h) - kiXu l,))\I,\ + — 

an ^-^ an 
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and refining the Riemann sums, 

— / {HX2, loig)) - fc(Ai, lo{g))) dg - — 1/| < lim inf ^ 
ctTT J J air c— i-oo c 

< lim sup ^ < — / (^(A2, /o(^)) - A;(Ai, lo{g))) dg + 

c^oo c air J J air 

The assertion of Theorem 1 foUows, as e > is arbitrary, and k{X2, lo{g)) = k{Xi, loiQ)) 
{g € (0, 00) \ /). 

This concludes the proof of Theorem 1. 
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